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ABSTRACT 

A great deal of evidence has recently been gathered in favor of the picture that Soft Gamma Repeaters and 
Anomalous X-Ray Pulsars are powered by ultra-strong magnetic fields (B > 10 14 G; i.e. magnetars). Nevertheless, 
present determination of the magnetic field in such magnetar candidates has been indirect and model dependent. 
A key prediction concerning magnetars is the detection of ion cyclotron resonance features, which would offer a 
decisive diagnostic of the field strength. Here we present the detection of a 5 keV absorption feature in a variety 
of bursts from the Soft Gamma Repeater SGR 1806-20, confirming our initial discovery (Ibrahim et al. 2002) and 
establishing the presence of the feature in the source's burst spectra. The line feature is well explained as proton 
cyclotron resonance in an ultra-strong magnetic field, offering a direct measurement of SGR 1806-20's magnetic 
field (B w 10 15 G) and a clear evidence of a magnetar. Together with the source's spin-down rate, the feature also 
provides the first measurement of the gravitational redshift, mass and radius of a magnetar. 

Subject headings: Pulsar: Individual (SGR 1806-20) — Stars: Magnetic Fields — Stars: Neutron — Stars: 
Magnetar — X-Rays: Bursts — Gamma Rays: Bursts 



1. INTRODUCTION 

Soft Gamma Repeaters (SGRs) are a rare, enigmatic type of 
pulsars that are distinguished by the hallmark of episodic emis- 
sion of short (~ 0.1 s), hyper-luminous (~ 10 39 - 10 42 ergs" 1 ) 
bursts of soft 7-rays (see e.g. Hurley 2000 for reviews). Less 
frequently, SGRs also emit giant flares that last for hundreds 
of seconds, exhibit remarkable pulsations, and even have mea- 
surable effects on the earth's outer atmosphere (Mazets et al. 
1979; Hurley et al. 1999; Ibrahim et al. 2001; Inan et al. 1999). 
The high flux, short burst durations and lack of evidence for ei- 
ther a binary companion or an accretion disk inspired appeal to 
unusual energy sources for SGRs, including objects with ultra- 
strong magnetic fields (Duncan & Thompson 1992), accretion 
onto neutron stars from a fossil disk (Marsden et al. 2001), and 
emission from strange quark stars (Zhang et al. 2000). 

SGR 1806-20 underwent two major bursting episodes since 
its discovery in 1979, the latest of which occurred in 1996, 
when Rossi X-ray Timing Explorer (RXTE) observations re- 
vealed a 7.47 s pulsar associated with the source that spins 
down at a relatively fast rate of ~ 10" 11 s/s (Kouveliotou et 
al. 1998). The detection of pulsations and high spin-down 
rates in the quiescent emission of SGR 1806-20 and SGR 
1900+14 (Kouveliotou et al. 1998; 1999), when interpreted 
within the magnetic dipole model, gave strong indication that 
these sources are endowed with ultra-strong magnetic fields 
(B > 10 14 G) for which the name magnetars was coined. How- 
ever, the lack of independent spectral confirmation of the mag- 
netic field strength and the irregularities in the spin-down rates 
allowed for alternative models (Marsden et al. 1999). 

Recently, a 5 keV absorption feature was detected in a 
burst from SGR 1806-20, providing evidence for proton cy- 
clotron resonance in an ultra-strong magnetic field (B w 10 15 
G) (Ibrahim et al. 2002). Here we present a comprehensive 
study of a large number of bursts from the source observed with 
RXTE during its last active phase in November 1996. We detect 
the 5 keV feature with high confidence in the overall data set. 
In light of the feature's characteristics and the source's obser- 
vational properties, proton cyclotron resonance is the plausible 



interpretation, providing the first confirmed evidence of pro- 
ton/ion spectral signature from a cosmic source. We use the 
feature together with the spin-down rate to further constrain the 
magnetic field and to obtain estimates on the gravitational red- 
shift, mass and radius of the SGR. 

2. OBSERVATION AND DATA ANALYSIS 

We obtained the RXTE Proportional Counter Array (PCA) 
data from the HEAS ARC archives and used the event-mode and 
good-xenon data to construct the burst light curves and spec- 
tra. We detected over 100 bursts with a variety of strengths that 
range from faint bursts, which do not have enough counts to ac- 
cumulate a meaningful spectrum, to very bright bursts, which 
have large dead-time and pile-up effects and in some cases sat- 
urated the PCA. For cogent spectral analysis, we considered 
moderately bright bursts which gave enough counts for accu- 
rate spectral fitting but at the same time have count rate not 
exceeding 90,000 counts/s so that dead-time and pile-up effects 
would not be significant (Ibrahim et al. 2001). Although such 
instrumental effects are not known to produce artificial spectral 
features but would rather harden the overall spectral continuum 
and modify the widths and depths of real features, this careful 
selection ensures results that are intrinsic to the source. 

We classified the bursts according to their temporal profiles 
into single-peak and multi-peak bursts and fitted each emission 
peak separately. This resulted in a total of 56 spectra 1 . For 
background subtraction, we used quiescent pre-burst data seg- 
ments. We found qualitative characterization of the spectral 
continua with an absorbed power-law model, which provided 
slightly better fits than an absorbed bremsstrahlung model for 
the overall data set. The fits showed a number of spectra with 
relatively large y 2 and narrow absorption dips in the vicinity of 
5 keV. To assess the significance of the dip, we added Gaussian 
absorption lines to the fit models and evaluated the change in 
X 2 using the F-test. The fits improved considerably, giving rise 
to a highly significant (> 4er) 5 keV absorption feature in some 

1 All burst observations have the same pointing coordinates and PCA EDS 
configuration and were analyzed identically. 
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bursts while others showed only modest improvement (< 3 a). 
Replacing the Gaussian line with an absorption edge was less 
significant by an order of magnitude. 

Fig. 1 shows examples of the bursts with strong evidence of 
the feature. The spectral characteristics of the burst continua 
and the line features are shown in Table 1. All features are 
narrow; 90% confidence upper limit on the line widths range 
from 0.45 keV to 0.95 keV. The features are fairly strong with 
equivalent widths (EW) varying from 400 eV to 850 eV. 

In Fig. 2 we show one of the bursts that did not show the 
feature significantly, for which we place a 90% confidence up- 
per limit on the equivalent width of 200 eV. The fact that these 
bursts have similar temporal profiles and count rates compara- 
ble to those in Fig. 1, yet show no features argues against in- 
strumental effects. In fact, the equivalent widths of the features 
shown in Fig. 1 exceed by an order of magnitude those of pos- 
sible artificial features that could be due to imperfections in the 
detector response matrix. As an additional verification, we ana- 
lyzed Crab pulsar spectra of similar counts using an observation 
that took place just a few hours before those of SGR 1806-20. 
The spectrum is well represented by an absorbed power-law and 
the fit is not improved by adding an absorption line near 5 keV. 
The 90% upper limit on the equivalent width is 40 eV. 

No evidence is found for a correlation between the appear- 
ance of the feature and spectral hardness, burst duration, counts, 
or peak count rate. Interestingly, most of the spectra that 
showed the feature significantly were associated with average 
2-50 keV flux in a rather narrow range of (6.34-9.66) x 10~ 8 
erg cm -2 s" 1 (the flux varied among all spectra between (4- 
200) x lO^ergcm" 2 s" 1 ). 

To estimate the overall significance of the 5 keV feature in 
the entire data set, we individually fitted each burst spectrum 
with an absorbed power-law and obtained the cumulative \\ 
and degrees of freedom v\ . Next, we added a narrow Gaussian 
line whose energy is allowed to vary around 5 keV only within 

1 keV, as implied by the 90% confidence level on the features 
in Table 1, and fitted the data again and obtained the new cumu- 
lative xl an d V2. We then performed the F-test on the two x 2 
distributions. With (x\ = 2630.3, v x = 2176) and (xl = 2302.3, 
h>2 = 2063), we get an overall chance probability p = 2 x 10~ 16 . 
If we exclude the five bursts shown in Fig. 1, we still detect the 
feature significantly with a chance probability of 2 x 10" 5 . 

Alternatively, if we assess the overall significance of the fea- 
ture from the 5 bursts in Fig. 1 using 2 p = n , ^[ nV pi pi P3 ■ ■ ■ p n (1 - 

Pmin) {N ~ n \ we get p = 5 x 10" 14 . This is a conservative estimate, 
not accounting for the presence of the feature in other bursts. 

3. DISCUSSION 

With the growing evidence of ultra-strong magnetic fields 
in SGRs (Kouveliotou et al. 1998; 1999) and AXPs (Hulle- 
man et al. 2000; Kern & Martin 2002), several models have 
recently been developed for proton and ion cyclotron features 
(Bezchastnov et al. 1996; Zane et al. 2001; Ho & Lai 2001; 
Ozel 2002; Thompson 2002a). Recent high-resolution obser- 
vations of SGRs and AXPs quiescent emission revealed no ev- 
idence for spectral features (Kouveliotou et al. 2001; Fox et al. 
2001; Patel et al. 2001; Juett et al. 2002; Tiengo et al. 2002). 
However, two spectral lines at 6.4 keV and 5.0 keV have re- 
cently been detected in two bursts from SGR 1900+14 and 
SGR 1806-20, respectively, suggesting that spectral features 

2 N is the total number of spectra inspected (56), n = 5, Pi-.-Ps are the chance 
probabilities given in Table 1, and />„„■„. is the least significant p. 



could be a property of the bursts alone (Strohmayer & Ibrahim 
2000; Ibrahim et al. 2002). While the 6.4 keV feature did not 
provide conclusive evidence on the magnetic field strength of 
SGR 1900+14, the 5.0 keV feature gave evidence for proton 
cyclotron resonance with an implied magnetic field of w 10 15 
G for SGR 1806-20. Our new results reported here confirm the 

5 keV feature in SGR 1806-20 burst spectra. 

The possibility that the 5 keV feature is due to electron 
cyclotron resonance in a typical neutron star magnetic field 
(B « 10 12 G) is discussed in Ibrahim et al. 2002. That the line 
width upper limits are up to an order of magnitude smaller than 
both observed and expected widths of electron cyclotron fea- 
tures argues against this interpretation. Moreover, an electron 
resonance at 5 keV would give a field of ~ 6 x 10 11 G, three or- 
ders of magnitude less than that derived from the spin-down rate 
estimate (~8x 10 14 G). The possibility that the line could be 
formed at several neutron star radii from the surface is unlikely 
due to the presence of a range of field strengths that would not 
generate a narrow feature. The splitting of Landau levels in 
B~6x 10" G (Pavlov et al. 1991), which seems potentially 
promising to distinguish an electron feature from an ion signa- 
ture if observed with a high resolution spectrometer, might be 
overwhelmed by Doppler broadening. Actual measurement of 
a narrow width would argue for the ion cyclotron interpretation. 

On the other hand, the observed line width, energy, and 
equivalent width are consistent with model predictions for pro- 
ton cyclotron resonance proposed by Zane et al. 2001 and Ho 

6 Lai 2001 for fully ionized Hydrogen atmospheres. The mag- 
netic field implied from the proton interpretation (~ 1 x 10 15 
G) is in good agreement with the independent estimate (~ 
8 x 10 14 G) obtained from the spin-down rate during the same 
observations (Kouveliotou et al. 1998). Lai & Ho (2002) and 
Ozel (2002) pointed out that vacuum polarization effects would 
reduce the equivalent width of the proton feature in quiescent 
emission. Bulik and Miller (1997) studied the effect of vacuum 
polarization on the burst spectra, but not in connection with the 
proton resonance. Appropriate models for burst spectra that 
include both effects have not been worked out as yet. In con- 
clusion, on observational grounds, a cyclotron resonance at 5 
keV together with the observed period and spin-down rate of 
the neutron-star source are consistent only if the star's magnetic 
field is greater than 10 14 G. 

In addition to the proton resonance at E p in B > 10 14 G, 
a-particles and ions with Z/A =1/2 would resonate in the 
same field at E p /2. However, the significant increase in the 
binding energy of atoms in a strong magnetic field (E (Z) oc 
Z 3 [In (B/Z 3 ) + const.] 2 , Duncan 1999; Lai 2001) would likely 
suppress the ion features. If the 5 keV feature arises from ion 
and a-particle resonance in B <~ 2 x 10 15 G, one would expect 
a similarly strong feature at <~ 10 keV due to proton fundamen- 
tal and ct-particle/ion first harmonics. This is not seen. There 
is a hint of a feature near 10 keV in some bursts, but it is more 
consistent with being a 5 keV second harmonic. Evidence of 
higher harmonics was seen in the first burst that showed the 
feature (Ibrahim et al. 2002). 

Using the proton fundamental resonance formula E p = 
6.3(l+z) _1 (B/10 15 G) keV 3 , and the observed resonance val- 
ues, the magnetic field in SGR 1806-20 is directly determined 
by B = (7.7-9.4)(l+z) x 10 14 G. For a given E p and B, the 
resonance relation gives a locus of M and R. Similarly, the 

3 Here z is the gravitational redshift and (1 + z)~' = (1 -2GM/Rc 2 )'/ 2 , and M 
and R are the mass and radius of the star, respectively. 
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Table 1 

Spectral Characteristics of some SGR 1806-20 Bursts with Strong Evidence For the 5 keV Absorption Features 



Burst 


Fluence 
(counts) 


Photon Index 


(10 22 cnT 2 ) 




E 

(keV) 


EW 

(eV) 


Width 

(kev) 


Xl/V2 


Pp-test 


Burst 1 


301 


1.43 ±0.25 


20.0 ±7.95 


48.2/32 


5 35+0 21 


850 ±40 


< 0.47 


25.3/30 


6.3 x 10~ 5 


Burst 2 


1422 


2.11 ±0.1 


22.7 ±2.53 


80.2/43 


3zy -0.15 


560 ± 110 


< 0.75 


44.3/41 


5.2 x 10"" 


Burst 3 


1891 


1.68 ±0.08 


20.3 ±2.19 


58.1/43 


5 69 +lul 


430 ± 100 


< 0.8 


37.7/41 


1.4 x 10-" 


Burst 4 


625 


1.47 ±0.14 


18.5 ±4.12 


57.6/43 




650 ± 170 


< 0.45 


37.4/41 


1.4 x 10~ 4 


Burst 5 


1914 


1.77 ±0.08 


20.4 ±2. 15 


66.5/43 


5 16^ us 


400 ± 100 


< 0.95 


49.6/41 


2.5 x 10" 3 



NOTE. — Errors quoted are the lcr error, except for the energy where the 90% confidence level is given. The width values are the 
90% confidence upper limit. Burst 1 is the shortest burst detected with the smallest number of background subtracted counts. With the 
C-statistic (Cash 1979), the random probability is 2.7 X 10~ 3 . 




Energy (keV) Energy (keV) Energy (keV) 

FIG. 1 . — Bursts from SGR 1806-20 with strong evidence for the 5 keV spectral feature. The first column panels show the burst light curve profile. The second 
show the fit residuals with an absorbed power-law model. The third and fourth show the spectrum and residuals, respectively, of the best-fit model including a 
narrow Gaussian absorption line. The residuals are shown in units of <r. The 5 keV absorption feature is present in the whole of bursts 1 and 3 and in part of bursts 
2, 4, and 5. The bursts times are 50405.3704, 50405.5051, 50405.3334, 50405.1843, 50405.4675 (MJD), respectively. 




Energy (keV) Energy (keV) 



FIG. 2. — An example burst with modest or no evidence for the 5 keV feature. The two-peak burst shows random residuals in the first peak and weak absorption 
excess (2.4cr) in the second peak. The burst time is 50405.5191 (MJD). 
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magnetic field-spin-down relation 4 , B 2 = y^t^z jp, gives an- 
other locus of M and R for a particular B. The intersection of 
the two relations constrains the mass and radius of the SGR 
as illustrated in Fig. 3. While the dipole model comes from 
an idealized case, it provides adequate order of magnitude es- 
timate for the magnetic field. Also, since the connection of 
SGR 1806-20 and SNR G10.0-0.3 is now uncertain (Gaensler 
et al. 2001), there is no longer a requirement that it has a large 
wind of particles to fill a plerionic nebula. It is then possible 
for the magnetic dipole radiation to be the dominant source of 
spin-down. 




0.5 

o.d F - 

8 10 12 14 16 

R (km) 

FIG. 3. — Mass-Radius curve for SGR 1806-20 as obtained from the in- 
tersection of the proton cyclotron resonance relation (dashed lines) and the 
spin-down rate relation (dotted curves) for an oblique rotor model. Both rela- 
tions are plotted at B = (1.0— 1.3) X 10 15 G (increasing from bottom to top). 
The line along the intersection points corresponds to the best-fit cyclotron res- 
onance energy while the other two refer to the 90% confidence level. Here we 
used the feature with the largest redshift (E p = 5.16 keV) and assume that it 
comes from very close to the stellar surface. For equations of state of current 
interest, £ would be in the range 0.3-0.45 for 0.1 < (M/M )/(i?/km) < 0.23 
(Lattimer & Prakash 2001). 

Assuming an oblique rotator and the average value (0.5) of 
sin 2 #, the M-R constraints above, along with the equations 
of state, would limit the radius and mass of SGR 1806-20 to 
R = (10.5 - 12.3) km and M = (1 .3 - 1 .8) Af if the neutron star 
were composed of nucleons and hyperons. The corresponding 
gravitational redshift and mass-radius ratio are z = 0.26-0.33 
and M/R = (O.12-O.15)M /km, which yield a surface mag- 
netic field B = (0.97-1.3) x 10 15 G. We note that our estimates 
of z and M/R are close to those recently reported from an iso- 
lated neutron star that showed 0.7 and 1 .4 keV absorption fea- 
tures (Sanwal et al. 2002). For sin 2 f? = 1, the M-R locus moves 
to values of R < 10 km, too small for equations of state that do 
not have quarks or other 'exotic' particles. A quark star would 
give R = (8 - 10) km and M = (1 .3 - 1 .95) M . 

With our present knowledge, there are several mechanisms 
that could explain the observed features in light of the magnetar 
model (Thompson & Duncan 1995). SGR bursts have been 
identified with photon-pair plasma that is created and trapped 
in the magnetosphere as a result of magnetically driven crust- 
quake. The plasma dissipates its energy by X-ray emission; part 
is radiated away as the burst and part is radiated back to the 
exposed surface of the neutron star. In this framework, proton 
cyclotron resonance could arise in several scenarios. 

4 P and P are the period and spin-down rate, 6 is the angle between magnetic 
dipole and rotation axis, and £ is the neutron star's moment of inertia in units 
of MR 2 . 



First, it could originate near the surface of the neutron star, 
in the local region of the crust-quake, where hydrogen in the at- 
mosphere is being irradiated by the plasma photons. Secondly, 
the cooling phase of the burst emission is expected to involve 
an ablation process that blows off baryons and ions from the 
surface of the star to its atmosphere. If this occurs in the first 
peak of a multi-peak burst (or in a single-peak burst), the opti- 
cal depth during a subsequent peak (or a later burst) is likely to 
be larger, thus augmenting resonance features. 

Features produced by these mechanisms come from near the 
surface and thus can be obscured by the rotational orientation 
of the star if the emission region is not along our line of sight. 
Note that the durations of the absorption line episodes are much 
smaller (< 0.01 P) than the pulsar rotational period. It was re- 
cently pointed out that emission from near the equatorial plane 
will be pronounced at the cyclotron resonance, in contrast with 
that from the poles (Thompson et al. 2002b). 

A third origin of proton resonance is from the baryon-loaded 
sheath, which constitutes the plasma outer shell. The resonance 
could occur in that region as protons in the baryon sheath are 
illuminated from within by the hard burst photons. 

In all the scenarios above the burst is essential, giving rise 
to optical depths and fluxes of high energy photons far exceed- 
ing those in quiescence, which may explain the absence of the 
feature from SGR and AXP persistent emissions. The effect of 
vacuum polarization, together with the correlation with the flux, 
the rotational and latitude constraints mentioned above, and the 
fact that the optical depth to resonant cyclotron scattering is 
proportional to the local density of protons could in principle 
account for the transient nature of the feature and its absence or 
weak presence in some bursts. 

In summary, the discovery of the 5 keV absorption feature 
provides a unique opportunity not only to diagnose SGR 1 806- 
20 magnetic field but also to reveal insight on the gravitational 
field strength, mass and radius of SGRs. This is an important 
step toward understanding a mysterious constituent of the Uni- 
verse. Further observations and careful analysis could allow 
unprecedented studies of physical phenomena in the extreme 
limits of gravitation and magnetic field, uniquely found in mag- 
netars. 
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NAG5-1 1951 and The George Washington University. 
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